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ABSTRACT: In this paper we use one-loop chiral perturbation theory in order to com- 



This makes it possible to investigate three systematic effects that plague lattice computa- 
tions: quenching, finite- volume effects, and the fact that lattice computations have been done 
at unphysical values of the quark masses and pion external momenta (only this latter effect 
shows up at tree level). We apply our results to the most recent lattice computation (ref. 
[1]), and find that all three effects are substantial. We conclude that one-loop corrections 
in chiral perturbation theory help in explaining the discrepancy between lattice results and 
the real-world value. We also revisit Bj{, which is closely related to the K + — > n + n decay 
amplitude by chiral symmetry. 
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1. Introduction 

One of the aims of practitioners of lattice QCD is the accurate computation of QCD weak 
matrix elements. There are two reasons for this. First, such matrix elements are needed in 
order to explain the experimental values of certain quantities. A prominent example is the 
A/ =1/2 rule for nonleptonic kaon decay. Second, precise knowledge of the values of weak 
matrix elements helps with the determination of less well-known entries in the Cabibbo- 
Kobayashi-Maskawa matrix. An important example is the kaon B-parameter B^, which 
parametrizes the strength of CP-violation in K — K mixing. (For these and various other 
phenomenological aspects of kaon weak interactions, see ref. [2] and references therein; for 
the lattice approach, see refs. [3,4,5] and references therein.) 

Actual lattice calculations [6,7,8,9,1] are hampered by many systematic errors. Even if 
the continuum limit is obtained, several important sources of error still remain. Two of these 
are caused by the use of a finite volume and by the use of the quenched approximation. A 
third systematic effect typically occurs for matrix elements involving more than two external 
particles. Taking K — ► 2n decay as an example, the matrix element most easily determined 
from a numerical computation is that with degenerate quark masses, and with all external 
particles (both the kaon and the two pions in this case) at rest. This does not correspond to 
the physical situation, where of course mj{ > 2m 7r , and each of the final pions has a spatial 



momentum |k| = \ \rn\ — m% (in the kaon's rest frame). This implies that such lattice 
results will have to be extrapolated to physical values of the quark masses and the external 
momenta. 

All of the above described effects can be studied in chiral perturbation theory (ChPT). 
Finite- volume effects can be calculated in ChPT following ref. [10], while a quenched for- 
mulation of ChPT ("quenched chiral perturbation theory," or QChPT) was developed in 
ref. [11,12]. The third systematic effect introduced above can also be investigated in ChPT 
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by simply calculating both the "unphysical" matrix element (with all external particles at 
rest) and the "physical" one (with the correct kinematics). The ratio of the two then gives 
an estimate of the factor needed in order to convert the unphysical matrix element into the 
physical one. This was done at tree level for the decay K + — ► n + ir in refs. [9,1] (at tree level 
no finite- volume or quenching effects show up in ChPT). It was found that even with the 
conversion factor (which also corrects for the fact that the quark masses in numerical com- 
putations are different from the physical quark masses), the lattice value of the K + — ► n + n 
matrix element is about two times larger than its experimental value. 

In this paper we study the K + — ► n + n weak matrix element and B^ in ChPT to one- 
loop order. This is interesting for various reasons. Nontrivial finite-volume and quenching 
effects show up in ChPT first at one loop, and therefore a one-loop calculation will give us 
insight into the size of these effects. Also, there are one-loop corrections to the tree-level 
conversion factor between unphysical and physical K + — ► n + n matrix elements, and one 
would like to know whether these corrections can be part of the explanation that numerical 
results appear to come out too large [4]. Moreover, it has been suggested that one can expect 
such one-loop effects to be sizeable, since they include an estimate of the effect of final-state 
interactions of the two pions [13]. 

We restrict ourselves to Bj^ and K + — ► n + n decay because these two quantities require 
the introduction of only one new parameter at the lowest order in ChPT. This is because the 
corresponding weak operators are both components of the same 5'?7(3)-flavor 27-plet, and this 
27-plet has a unique representation in lowest order ChPT. It would of course be interesting to 
also consider K — ► 7r + 7r~ decay, but this would lead to the introduction of new parameters 
(corresponding to the SU(3) octet operators that contribute to this decay). Also, it has been 
observed that inclusion of the a resonance may be required for a full understanding of this 
decay on the lattice [14], which would take us outside of the systematic approach of ChPT. 

Beyond lowest order in ChPT, a large number of new operators exist which can contribute 
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to B^ and to K + — ► tt+tt decay [15], each with its accompanying new parameter (the so- 
called 0(p ) low-energy constants). The values of these new parameters are largely unknown. 
Therefore, at one loop we can only estimate the size of the nonanalytic terms (using a 
reasonable value of the cutoff), because the nonanalytic terms do not depend on these new 
parameters. In this paper, we will assume that using a reasonable value of the cutoff in 
the nonanalytic terms and setting the 0(p ) low-energy constants to zero will lead to a 
valid estimate of the size of corrections to tree-level ChPT results. Some idea about the 
quantitative uncertainties introduced by this approach can be obtained from the sensitivity 
of our results to the value of the cutoff. These issues will all be discussed in more detail in 
due course in this paper. 

Let us end the introduction with an overview of the paper. We revisit B^ in section 
3 (quenching and finite-volume effects for B^ have been calculated within ChPT before in 
refs. [12,5]), while the calculations for the K + — ► n + ir matrix element are contained in 
sections 4-6. In section 3 we also discuss our approach to the use of ChPT at one loop 
in some more detail. In section 4 we present our calculation of the nonanalytic one-loop 
corrections for the physical matrix element. Section 5 explains which Euclidean correlation 
function is computed in numerical lattice QCD, and in section 6 we present in some detail 
the calculation of this correlation function in one- loop (Q)ChPT. For this calculation we 
restrict ourselves to the case of degenerate quark masses (m u = m^ = m s ). In section 7, 
we use the results of the previous sections in order to obtain quantitative estimates of the 
various systematic effects for typical lattice computations. We start with calculating the size 
of the nonanalytic one-loop corrections in the real world, as was done before by Bijnens et 
al. in ref. [16] and by various other authors [17,18,19,20]. We then estimate the one-loop 
correction to the conversion of the unphysical lattice K + — ► n + n matrix element to the 
physical one in the unquenched theory. Next, we include the changes that result from using 
the quenched approximation for the unphysical matrix element, and we discuss the numerical 
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results of ref. [1]. The uncertainties introduced by our lack of knowledge of all low-energy 
constants are discussed. We end this section by considering a quantity which is basically 
the ratio of Br and the K + — ► n + n matrix element [21]. This quantity is interesting 
because it is independent of the tree-level low-energy constant, and therefore not sensitive 
to ambiguities in its value. All calculations are repeated for the finite spatial volume which 
was used in ref. [1]. Our conclusions are in section 8. Section 2 contains a summary of 
other (Q)ChPT results that we will need in this paper. Appendix A presents an adaptation 
of the quenched Br result relevant to staggered fermion QCD, and Appendix B relates our 
result for the K + — ► n + n correlation function to the general analysis of such Euclidean 
correlation functions in ref. [22]. Appendix C deals with a technicality. 

2. Essentials of ChPT and QChPT 

We start with a summary of results in ChPT [23,24,25] that we will need in this paper. 
This also gives us an opportunity to establish our notation. We will refer to the theory 
without quenching as the "full" or "unquenched" theory. The SU (3)-octet Goldstone meson 
fields are organized in the nonlinear representation 

£ = exp — , (1) 
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where <b is the hermitian matrix 
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and / is a parameter with the dimension of a mass. Masses of the u, d and s quarks are 
incorporated in the mass matrix 

m u o o 

M = | m d j (3) 

m s 



To lowest order (which we will refer to as 0(p )) in ChPT, the Euclidean Lagrangian is 
constructed with the introduction of an additional parameter v : 

c2 



£ = i-Trfd^Edfj:^ - vTi(mt) + M+£ 



(4) 



One can read off the tree-level meson masses from the quadratic terms, after expanding E 
in (for simplicity we will choose m u = m^ = m): 
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It follows that 
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which predicts a value of the r/-mass about 3% too large. 

At tree level, the parameter / is the weak decay constant, / = f n = fj{. At one loop the 
pion and kaon decay constants f n and fg and corresponding wavefunction renormalizations 
Z n and Zj( are given by 
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for the kaon, where 
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In these equations, we made use of the fact that we can replace the tree-level quantities rrv^ q, 
mj^Q and tttc g by their renormalized physical values m^, m^ and mi. The parameter y n 
is defined as 

V* = ^ ■ ( 13 ) 

A is the cutoff, and the dimensionless parameters K\ and K<± are contact-term coefficients 
coming from the 0(p ) chiral Lagrangian [25]. 

A Lagrangian formulation of quenched ChPT (QChPT) was developed in ref. [11]. The 
basic feature is the introduction of (unphysical) spin 1/2 "ghost" quarks with bosonic statis- 
tics and the same masses as the physical quarks, so that their virtual quark-loop contributions 
cancel those from the original quarks [26]. As a result, the 77 meson is light [27,11,12], and 
therefore needs to be kept in the chiral Lagrangian. The original octet Eq. (1) has to be 
extended into a nonet to include the 77 : 



£ g = e^/ V3/ <>£, (14) 

where f q is the quenched parameter equivalent to / in the full theory, and / is replaced by fq 
in Eq. (1). In the mass-degenerate case m u = m^ = m s , the 77 two-point function acquires 
both single- and double-pole terms: 

(9 9 \ 

jrr^-- ,2 + T,2 J ; ( 15 ) 

where fi is a parameter which would correspond to the singlet part of the 7/ mass in full QCD. 
In full QCD, an estimate using the experimental mass of the 77' gives: /j, /3 ~ (500 MeV) 
(for a — 0). In the mass-nondegenerate case m u = va^ 7^ m s , the 77 two-point function also 
acquires a double-pole term through mixing with the 77 . It is then simpler to work in the 
(nondiagonal) basis of uu, dd, and ss meson states in the neutral sector (labeled by i — u, d, s 
respectively). Their two-point functions are: 

G,JP, ,) = S(P + ,) (-^ - (p2 f^+M^ I ' (16) 
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where M 2 = Svnrn/ f 2 When m u = m^ ^ m s , it is easy to see that M 2 
Mg = 2m 2 K - m 2 Q . We define 
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for later use. 



The quenched one-loop expressions for the weak decay constants f n and fj{ and wave- 
function renormalizations Z n and Zj{ in the mass-nondegenerate case (m u = m^ ^ m s ) 
are 



fn = fa 1 



Zjf 



m: 



(4vr/ g ) 2 



K 



2m, z 



3(4tt/ ? ) 2 



K 



fK = fq[l+lK 



ill 



K 



(W,) 2 



K 



iff 



3 



77? 



«jr 



if 



(4vt/,) 2 



K 



(18) 
(19) 
(20) 
(21) 



where 



-8[ 1 + — ^ : 

2 V 2(1 -y n ) 



log 



Vn 



JJtt 



III 



-a 



K 



3 (Aixf q 



, . 2/tt(2 — 2/tt) , 

1 + ^71 r log 

2(1 — 2/tt) 



JJn 



Vir 



(22) 



and the dimensionless parameter K is again a contact-term coefficient from the 0(p ) La- 
grangian. In contrast to the full theory, there is only one such coefficient in the quenched 
case. In ChPT, an 0(p ) operator containing the trace of the mass matrix M, Eq. (3), when 
evaluated at tree level, contributes to a particular linear combination of K\ and Ki in the 
unquenched quantities. However, in QChPT, the same operator is converted into one where 
the supertrace [11] of the mass matrix, now enlarged to include the ghost-quark masses, is 
taken, and hence its contribution to Eqs. (18,19,20,21) vanishes. Note that chiral logarithms 
and 5 or a dependence are absent in Eq. (18) and Eq. (19). This is a consequence of taking 
m u = m d . 



In general, the double poles in Eq. (16) lead to severe infra-red divergences in physical 
quantities in QChPT. In fact, physical quantities with a singular chiral limit and/or anoma- 
lous volume dependence were found in refs. [11,12,28,29,30,31,32]. An example is fx in 
Eq. (20), which diverges in the limit y^ — ► 0. 

3. Review of Br 

The kaon B-parameter Br is defined as 

(K°\(sdsd) LL \K°) 

BK = 87232 ' ( 23 ^ 

V K m K 

in which the four-quark operator is defined by 

(Q1Q2Q3Qa)lL = (Q1L1^Q2l)(Q3L1^<14l) , (24) 

where q^ = ^(1 —75)5 is a left-handed quark field. The denominator in Eq. (23) is the matrix 
element (K \(sdsd) ll\K ) evaluated by vacuum saturation. In ChPT, weak-interaction 
operators, the chiral transformation properties of which are dictated by the Standard Model, 
are constructed from the Goldstone meson field £ and the mass matrix M . The operator 
(sdsd)LL is a component of a 27-plet of SU(3)l and has AS = 2. To lowest order in ChPT 
(0(p )), there is only one 27-plet operator: 

O' = a 27 ^(W^MW^ji , (25) 

where the tensor t£, satisfies: 

E4 = ' E*« = ° and % = t ii- ( 26 ) 

i i 

In order to select the AS = 2 component, we set 

41 = 1 (27) 



and all other components to zero. The parameter a 2 7 is not constrained by any symmetry 
and its value is determined by QCD dynamics. At tree level (or order p ), one obtains 

(K°\O'\K ) = ^fml . (28) 

The factor m\ comes from a factor p\ where we take the external kaon momentum p K on 
shell, i.e. m^ here is the renormalized physical kaon mass. Therefore, 

B{ uU = ^ = B-f (29) 

(fx = f at tree level). 

In this paper, we are interested in corrections of order one higher in ChPT power counting 
[24,25]. We will refer to these as 0(p ) contributions. In general, they consist of one-loop 
corrections coming from the 0(p ) operators along with tree-level contributions from 0(p ) 
operators. A large number of 0(p ) operators with the symmetries of the operator in Eq. (23) 
can be constructed from X and M; see ref. [15] for a complete list of 0(p ) operators that 
are relevant to strangeness-changing nonleptonic weak-decay processes involving up to four 
Goldstone mesons. Each of them is associated with an arbitrary coefficient, which, just like 
a 2 7, cannot be determined within ChPT. Since at 0(p ) only tree-level contributions from 
these operators are needed, contact terms with arbitrary coefficients result. In contrast, one- 
loop contributions from 0(p ) operators generically give rise to nonanalytic functions of the 
quark masses with coefficients determined by the 0(p ) operators, and a cutoff A naturally 
arises. The dependence on the cutoff of the 0(p ) -operator coefficients can then be chosen 
such that physical quantities are independent of the cutoff. 

In order to determine the 0(p )-operator coefficients, many authors have used combina- 
tions of particle phenomenology and theoretical ideas such as large- N c arguments, resonance- 
dominance models, etc. (see for instance refs. [17,18,19,20,25,33,34]). A less phenomeno- 
logical approach is to find relations between different physical quantities that are indepen- 
dent of 0(p )- and 0(p ) -operator coefficients. An illustrative example is the combination 
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frjf-K I fJ of the pion, kaon and eta decay constants [25]; examples in QChPT can be 
found in refs. [28,35,32]. In this approach, physical predictions at order p are made without 
the model-dependent uncertainties mentioned above. These relations of course are precisely 
the Ward-identities associated with spontaneously broken chiral symmetry. Theoretically, 
an infinite number of such relations can be constructed. However, the large number of low- 
energy constants at order p makes this difficult in practice. This is certainly true when we 
only consider B^ and K + — >■ % + ir decay, both of which depend on a large number of 0(p ) 
low-energy constants. (A program to determine some of these coefficients from experiment in 
the m^ = approximation has been carried out in refs. [18,19].) In this paper, we will limit 
ourselves mainly to the nonanalytic one-loop corrections. With a reasonable choice of the 
cutoff, the magnitude of the corrections to the tree-level expressions for certain quantities 
can be estimated. 

Next, we present the results for the full and quenched K — K matrix elements of the 
operator O' to one loop. For the full case, we obtain 

(K°\o'\K y f = ^^K ^i + _^k ]^ I + Fl + F2 + F3 _ (F2 + 2Fs) yn + F3 y 2 



Uir + K2 )-f y * 



(30) 
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while for the quenched case 
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with 2/tt defined in Eq. (13) and S in Eq. (17). In QChPT, the operator (sdsd)^^ is represented 
by Eq. (25) with a 2 7 replaced by a new parameter a%. The mass relation Eq. (6) has been 
used to eliminate m^ in favor of m^ and mx, except in the r] chiral logarithm of Eq. (30). 
The dimensionless parameters F^ and Fj with i = 1,2,3 and j = 1,3 come from 0(p ) 
weak operators. They are parametrized in such a way as to obtain simple expressions for 
B^ below. There are just two instead of three independent contact-term coefficients in the 
quenched case: the reason is the same as that for the analogous difference between the full 
and quenched 0(p ) expressions for f n ^ and Z^ ^ in the previous section. These results do 
not include finite-volume corrections, which are exponentially small in m^L for large volume 
(where L is the linear dimension of the spatial volume) [12]. 

The results for B^ to one loop for respectively the full and quenched theories in the 
mass-nondegenerate case are 
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The parameter e is defined as 

2 
vat. m s — m , . 

e = 1 - -f = * = 1 - 3/tt - 38 

(In Appendix A we present a modified expression for 5^ which may be applicable when 
quenched staggered fermions are used at nonzero lattice spacing and flavor symmetry is not 
fully restored.) The parameter B* is defined in Eq. (29), and B q is defined from af 7 and f q 
by an equation similar to Eq. (29). As is clear from these expressions, the full and quenched 
theories are entirely different theories. In particular, there is no reason to believe that the 
coefficients of the 0(p ) and 0(p ) operators should be the same in ChPT and QChPT. The 
5 and a terms in Eq. (37) come from the double poles in Eq. (16). The full ChPT result 
Eq. (36) has previously been calculated in ref. [16] in the m^ = approximation and in refs. 
[36,19,29,20]; whereas the QChPT result Eq. (37) has previously been calculated in ref. [12] 
for the mass-nondegenerate case with a = 5 = 0, and in ref. [5] for the mass-nondegenerate 
case with a — 0. We note a sign difference in the 5 term between Eq. (37) and the result in 
ref. [5]. 

Note that, apart from the usual nonanalytic cutoff dependence, nonanalytic functions 
depending on e appear in Eq. (37). This dependence is uniquely determined by (Q)ChPT, 
and predicts the dependence of Bj{ on m^ and m^ once the values of the low-energy 
constants K^, K, F^ and Fj have been determined at a given value of the cutoff A from 
lattice results for B^ at given fixed values of these masses. Note that the contributions 
proportional to 5 and a in B*L vanish in the mass-degenerate case, e = 0. (This is why 
we have not absorbed contact terms proportional to a into F\ and F%.) We see that, apart 
from a change from B* to B q and / to f q , quenching does not introduce any change in the 
nonanalytic one-loop corrections of Bj^ in the mass-degenerate case. As we will see, such is 
not the case for the K + — ► n + n matrix element. 
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4. K + — ► 7r + 7r in the real world 

The AS = 1, A/ = 3/2 i^ + — > 7r + 7r decay amplitude is proportional to the weak 
matrix element 

(n tv | (sduu + smki — sdd<£) LL \K ) . (39) 

The operator is the Al = 3/2 component of the same 27-plet that also contains the operator 
(sdsd)^^ (Eq. (23)). To lowest order in ChPT, the A/ = 3/2 0(p ) -operator is represented 

by 

4 = a 27 rg(E^Et) ifc (E^Et) j7 , (40) 

where the tensor r, , has nonzero components 



21 12 12 21 1 

r 31 = r 13 = r 31 = r 13 = g ' 

HI = r§ = "^ (41) 

(all other components vanish). It is easy to check that the tensor r?, satisfies condition 
Eq. (26) . The parameter a 27 is the same as in Eq. (25). This was used in ref. [21] to 
derive a tree-level relation between the K + — »■ 7r + 7r° decay rate and Br: at tree level, using 
Eq. (29), we have 

{.+AOAK+) = ^|K " -4) = ^U (ml - ml) B K . (42) 

Again, the factor mj^—m^ arises from taking the kaon and pion momenta on shell, analogous 
to Eq. (28). 

Bijnens et al. [16] were the first to investigate Bj^ and K + — »■ 7r + 7r decay to one 
loop in ChPT. They used the approximation of vanishing pion mass based on the fact that, 
experimentally, ml/m K rs 1/13. Kambor et al. [18] gave a general analysis of strangeness- 
changing nonleptonic weak-decay processes involving up to four Goldstone mesons including 

K + — ► Tr + Tr , to one loop. 
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Figure 1. Feynman diagrams for the one-loop contributions to the K + — ► n + n matrix 
element of O4. Filled triangles denote weak-interaction vertices and filled squares denote 
strong-interaction vertices. 

We have repeated the calculation with nonzero pion mass and m u = m^ 7^ m s . Loop 
corrections come from diagrams (a) to (d) in fig. 1 and, in addition, from wavefunction 
renormalizations. Our result is: 
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The functions A(|/7r) and B{y^) are defined as follows: 

1 + y/l-lyn 



A(yn) 
B{yv) 



1 - 4y n log 



1 - Uy n + lQy% log 



1 - VI -%r, 

7 - 4^ + Vl- 44^ + 16^ 



(45) 
(46) 



7 - 4^ - Vl" 44^ + 16^ 
The parameter y^ is defined in Eq. (13). J z comes from wavefunction renormalizations 
(Eq. (8) and Eq. (10)), and It comes from converting / to f n in the tree-level contribution 
to Eq. (43), using Eq. (7). The associated 0(p ) contact terms with coefficients Kj have 
been combined. (Note that refs. [18,19] keep the bare parameter / in the tree-level term.) 
I a and Ifr are nonanalytic one-loop corrections from diagrams (a) and (b) respectively, while 
we have combined those from diagrams (c) and (d) in / c +^. The dimensionless parameters 
G\ and G^ come from 0(p ) weak operators [15]. We have again eliminated mi in favor 
of m^ and m K using Eq. (6) except in the r\ chiral logarithms. The imaginary term in 1^ 
constitutes the lowest nontrivial order phase information for K + — ► n + n decay from ChPT. 

Analytic continuation of both A(y n ) and B{y 1T ) is understood for values of y^ for which 
the arguments of the square roots would be negative. Of particular interest to us is the 
function B(y n ). The expression in Eq. (46) is applicable to the region extending from 
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y n — to the smaller of the two zeros of the argument of the square root, which is around 
0.023. In the region between this value and the larger zero, which is around 2.7, the following 
expression should be used: 



B(y w ) = -2yj-l + Uy n - IQyl arctan I v l ±^j- 1% ' j . (47) 

The real-world value for y n m 1/13 falls into this second region as opposed to the point 
y n — which belongs to the first region. 

If we do the calculation with m^ = m^, we find that the combined contribution of I a , 
1^ and I c +d vanishes. 

Setting m n = in Eq. (43), we get: 

(7r+7r |O 4 |i^ + > = 

Ylia^mL { m\ / 9 m\ 4 7 YiK\ Ko \ } , s 

^^f{ 1+ (S7p(-2 l0 «# + ^ + 3 ^3-2 + V-T + Gl )}( 48 ) 

The chiral logarithm of ref. [16] appears, in addition to nonzero constant terms. The 
numerical significance of this will be discussed in section 7. 

5. K + — > vt + vt matrix elements from the lattice 

K + — ► 7r + 7r decay has been studied on the lattice by two groups [8,9,1,4], both of which 
used degenerate quark masses and the quenched approximation. On the lattice, a kaon K + 
at rest is created at time t — 0, using an appropriate bilinear quark operator. At a much later 
time t\ a lattice version of the AS = 1 and AI = 3/2 weak operator in Eq. (39) annihilates 
the K + and creates its decay products that include tt + and n . Finally, at a time t<i much 
later than t\, by a judicious choice of annihilation operators, the two pions, both at rest 
(see below), are picked among the products and annihilated. While the equations in this 
section apply to both finite and infinite spatial volumes, we will focus on the torus which is 
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the case of interest to lattice computations. The whole process can be expressed in terms of 
the following correlation function: 

C(t 2 ,h) = (0|7r+(t 2 )7r (t2)O 4 (ti)K-(0)|0) (49) 

= V / i w i A 0\e Ht H + (0)7v (0)e- Ht2 \n}{n\e H ^O 4 (0)e- H ^\m}(m\K-(0)\0) 
^— ' nm mm 

n,m s ' ' ' ' ' 

(50) 

< 2 »*i»o (0|7r+(0)7r (0)|7r+7r )(7r+7r |O 4 (0)|^ + )(^ + |^(0)|0) -E 27r (t 2 - tl ) -m K t, 

{u+u q \u+it q ){K+\K+) 

(51) 



where H is the QCD Hamiltonian, vr + (t2) = fd x 7r + (x, t 2 ) and similarly for 7r (£2) an d 
If - (0). These fields couple to mesons with zero momentum. The operator O^iti) = 
fdx 4 (x, t\) couples only to states with zero total momentum. The state |7r + 7r ) is 
the zero-momentum two-pion state with lowest energy, and \K + ) is a zero-momentum kaon 
state. The sums over n and m run over complete sets of states. Eq. (51) shows the leading 
term of Eq. (50) in the large-time limits. Note that E2 n , which denotes the energy of two 
pions at rest in a finite volume, deviates from its infinite- volume value 277% [37,38]. 

We will also need the four-pion correlation function and the kaon propagator in order to 
extract the desired matrix element (n + n 10^(0)1 K + ): 

(0|7r+(t2)7r (t2)7r-(t 1 )7r (t 1 )|0) = ^ 7 l T (0|7r+(t2)7r (t 2 )|77)(n|7r-(t 1 )7r (t 1 )|0) 

z — ' (77 71) 

2 

(52) 



t^ m^jo^m^Ai 2 -EzAh-h) 



and 

(0\K+( tl )K-(0)\0} = J2j^-A0\K + (tl)\n}(n\K-(0)\0} 

z — ' (71 In) 

ti>0 \{K+\K- (0)\0)\ 2 u 

(K+\K+) ' [ ' 

Some comments are in order. First, after taking the square root of the factors in front of 
the exponentials in the large-time limits in Eq. (52) and Eq. (53), we can extract the desired 
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matrix element from Eq. (51) only if we know the normalizations of the state vectors. We 
choose the complete set of states to consist of single-particle eigenstates and their tensor 
products. We use the canonical relativistic normalization for single-particle states: 

(K+(p)\K+(q)) = 2L 3 E K (p)5 p , q , (54) 

where Ej^(p) is the energy of the kaon with spatial momentum p, L is the spatial volume, 
and <5 P; q is the Kronecker delta, with similar definitions for other mesons. The normalizations 
for two-particle states are 

(7r+(P)7r (p)|7r+(Q)7r°(q)> = (7T+(P)|7r+(Q)>(7r ( P )|7r°(q)> . (55) 

Second, the process corresponding to the large time leading order contribution in Eq. (51) 
does not conserve energy (except when mj{ = 2m i: ). We will therefore refer to the matrix 
element (ir + ir 10^(0)1 K + ) as "unphysical." The physical matrix element corresponds to an 
excited state in Eq. (50) if mj^ > 2m- K ^ sub leading in the large-time limits. In principle, 
since in this "unphysical" case O4 inserts momentum, its chiral representation could include 
total-derivative terms. However, there are no such terms to order p . It is clear that we are 
unable to uncover the phase of any matrix element, in accordance with the discussion in ref. 
[22]. 

6. Calculation of (tt + tt \O^(0)\K + ) in finite-volume (Q)ChPT 

While lattice computations are done at a nonzero lattice spacing (an issue we will not 
address in this paper), they are also confined to a finite four- volume. We will take the finite 
spatial volume as a cube with linear dimension L on each side. Gasser and Leutwyler, among 
others [10], argued that when periodic boundary conditions are used, spatial finite- volume 
effects can be taken into account in ChPT by changing the (spatial) momentum integrals in 
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loop computations into discrete sums, i.e. for any function /(k), the integral J ,i ^ /(k) is 
changed into the sum -X Sk/C 4 )' wnere k = 2im/L with n G -ZJ , the discrete momentum 
in the box. The prescription is valid up to corrections vanishing faster than any power of 
L . On the other hand, the values of the coefficients of all 0(p ) and 0(p ) operators 
remain the same as in infinite volume, again up to corrections vanishing faster than any 
power of L . 

A. Tree-level contribution 

We will now describe the calculation of C(t2,t\) (Eq. (49)), in finite-volume ChPT. We 
will choose t<i > t\ > 0. We first show the relatively straightforward tree-level calculation in 
order to set up notation that will facilitate the discussion of the one-loop calculation. We 
define the 7r + two-point function: 

D+(x,y) = (7r+(x,t)n-(y,s)) , (56) 

where x = (x, t) and y = (y, s) are space-time points. The two-point functions Dq and Dj^ 

for 7r and K + are defined in a similar way. Defining the sign function 

f 1, £>0 
sign(t) = <^ 0, t = , (57) 

[-1, £<0 

and 

S% = 1 , e+ = H+ = m* sign(t - s) , (58) 

we have for the spatial integrals of the two-point function and its derivatives: 

:0 



3,, n t„ „.\ - ^+ „— mJ\t— s\ 



d^D+(x,y) = -^e 
2m n 



jd^^D + (x,y) = 5„^±-e 



r 

— m-At— s\ 



(59) 



- : - e -m n \t-s\ 



with similar definitions of £q' s an d £,k s associated with the two-point functions Dq and Dj(, 
respectively. 
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At tree level, we substitute the relevant part of O^xli, ti) contributing to the K + decay: 



Aia 



4 (x 1 ,t 1 ) - -Jj^(4x-d> i iPd li K+ - d»ir-ir%K+ - 3^7r-^7r A'+)(x 1 ,t 1 ) , (60) 
where each field is taken at the same space-time point (x-i,ti), into Eq. (49) and get 



c(i 2 , ii) = --^ y ^x 2 y ^x 2 y ^ X1 y ^ 

AD + (x 2 , xi)8^Dq(x2, xi)d tl D K (xi,x K ) - d^D + (x 2 , xi)D (x 2 , xi)d IJi D K (xi,x K ) 

-3d^D + (x2,xi)d fi D (x 2 ,xi)D K (xi,x K )j , (61) 

where x\ = (xi,ti), x 2 = (x 2 ,£ 2 ), 5 2 = (x 2 ,£ 2 ), x K = (x#,0) and <9 M = g|^- . After 
integrating over spatial coordinates corresponding to the external particles, i.e. x 2 , x 2 
and x^, with the help of Eq. (58) and Eq. (59), an expression independent of xi is left. 
Integrating over x^ then gives a trivial volume factor L . The tree- level result is 

C(t2,h) = "^ ^ ^^{^+^K - &ZotK ~ 3£*5&J 

19?7v f 3 P -2m w (t2-*i) p-rriKh 

= 1^^ + ^L__£_. m 

Note that to this order i? 27r = 2??%. A comparison of mass prefactors in Eq. (62) and 
Eq. (42) already shows the difference between unphysical and physical matrix elements. We 
will extract the unphysical matrix element from C(t 2 ,ti) after the discussion of one-loop 
corrections. 

B. One-loop contributions in full ChPT 

In this subsection we discuss the one-loop contributions from diagrams (a) to (d) (see fig. 
1). We will from now on restrict ourselves to the mass-degenerate case. Diagram (a) is the 
tadpole diagram generated by only the weak operator 0^{t\). All octet mesons appear on 
the loop. Diagrams (b) to (d) involve the 0(p ) strong-interaction four-point vertex (£j n ^). 
These contributions are of the form 

C(i 2 ,£i)(b, c, d) = -(7T + (t 2 )7r (t 2 )(yrf 3 x s rft s £ mf (x s ,t s ))o 4 (t 1 )K-(0)) c , (63) 
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where (• • -) c denotes all connected Wick contractions, and where the relevant terms (cubic 
in the octet field (Eq. (2)) of the weak operator 0^(t\) are taken. They fall into different 
categories, represented by the different topologies of diagrams (b) to (d). Due to conservation 
of charge and isospin by the strong interaction vertex, only n + and n appear on the loop in 
diagram (b). In diagrams (c) and (d), a kaon and a pion or eta appear on the loop. We first 
explain in detail the contribution represented by diagram (b). The two external pions are 
Wick-contracted with two pions out of the four pion fields in Ci n p The other pion fields in 
dini are Wick-contracted with two pions from O4. The one field left in O4 is Wick-contracted 
with the external kaon. We will denote fields contracted between O4 and Ci n t by a subscript 
"/." For diagram (b), we end up with the product of 



(d^-^d v K + - d^njd^K+y^h) 



(64) 



and 



■ir-ir°d^irtduir T - d^'d^j^ + - 



l /^"J 



d v 



7T 7T 



7T 



-0M-O 



7T 



d^ + 7l°j + TT+fyTT? 



+ ml vr 7r°7r+7rJ Ux a ,t 8 j ■ 



(65) 



from O4 and Ci nt respectively, with a prefactor — Sia 2 T/(V2f ). Fields without subscript 
I contract with external lines. This product leads to a proliferation of terms with different 
possible combinations of derivatives acting on the fields, so let us represent each term in 
Eq. (64) and Eq. (65) symbolically by 



<r(xi,ti)<r|(xi,ti) and $ 2 (x s , t s )$|(x s , t s ) , 



(66) 



respectively. The fields a and $ stand for pion or kaon fields and their derivatives. With 
Eq. (66), all terms in Eq. (63) corresponding to diagram (b) can be written as 



A 2 d% Mq jd\ K d\ s fdt a (7T + (x 2 , t 2 )$(x s , t s )) (7T (x 2 , t 2 )$(x s , t 8 )) 



X 



($ 7 (x s , t s )o- 7 (xi, t 1 ))($ / (x s , t«)<T/(xi, t{)) x 
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dt s (fix s d\i 



(<r(xi,ti)K + (x K ,0)) 

e -2m 1T \t 2 -t s \ e -m K ti 



Am 2 2mx 

($ / (x, s ,t s ) ( x / (x 1 ,t 1 ))($ / (x s ,t s ) ( 7 / (x 1 ,t 1 )) . (67) 



The ^-factors are chosen according to the prescription of Eq. (59), and possible associated 
factors 8uA have not been displayed explicitly. At this point two simplifications occur. First, 
Lorentz contractions of derivatives on the same internal leg do not occur because of the fact 
that they come from the simple product of two Lorentz invariants O4 and L{ n t (see Eq. (64) 
and Eq. (65)). Second, if a derivative is Lorentz-contracted with another one acting on a 
field contracted with an external field, only the 4-th component of the derivative survives, 
according to Eq. (59). Using the momentum representation of the pion two-point functions, 
we obtain: 

f o p -2m n \t 2 -t s \ p -m K ti 



dt K I d X s / GTXl k x 

4m~. 2m k 



e -i(k-(x s -xi)+fc 4 (£ s -ii)) e -i(l-(x 8 -xi)+/ 4 (t a -ti)) 

k Jl (k 2 + m 2 ) (I 2 + m 2 ) 




where 



na 



M = —prk (-mtfsignfa)^ + k%) (k% - m^sign(t 2 - t s )i(k 4 + / 4 )) . (69) 

We use the abbreviation J*. = -p% J^ f -^ with k = 27in/L and n E 2Z . We then integrate 
over xi. This results in a delta function setting k = — l. Physically, this corresponds to the 
fact that 04(^1) only couples to states with zero total momentum. One of the momentum 
sums, say over 1, can be performed which leads to an integrand independent of x s , and 
subsequent integration over x s then gives a trivial volume factor L . The result is 

3 r e -2m v \t2-t 8 \ e -m K ti x r dkA r^ g-ik^-h) e -ih(t.-ti) 

L J dts 4^1 2m~^~ 1M-.J ^J2^ M (jfej + k 2 + mj) {I 2 + k2 + m 2 ) ' (?0) 
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where M now takes on the form: 



M 



na 



27 



V2f 



-m^sign(ii)i/c4 + /C4/4 — k j (^4/4 — k — m 7r sign(t2 — t s )i{k^ + Z4) J . 



The following formulas are needed for the evaluation of the integrals over k^ and Z4 in 
Eq. (70): 

dA; 4 fdk e -i(k4+U) (t s -h) 



where 



27T 7 27T (^+a#(^+U# 

' / ff 6(t 8 - h) - co k 5(t s - tl ) + ^ e -2«kl*.-*il , .A = fcjjj 
gn(t s - tO (>(t s - «0 - JJ e -2«kl*.-*il 



S1E 



v 4w k 



l_5(t s -tO-Je- 2 ^l*.-*il, 

_l p -2wk|t 8 -ti| 

^sign(t s -t0e _2a;k l* s_il l , 
l_ e -2w k |i s -ii| 



^k = V k 2 + m| . 



„4 = jfc|/ 4 



.A = A;| 

.4. = /C4/4 
A = ik^ 
A=l 



(71) 



(72) 



The strong interaction can "occur" at any time t s , a fact reflected by the integration 
over t s from —00 to +00. Carrying out the integrating over t s after integrating over k^ and 
£4, we finally obtain for the contribution from diagram (b): 



24ia 97 L 3 o e - 2m *(*2-ti) e - " 1 ** 1 



>/2/ 



27-^ 2 

— m 7T 



4m 2 



2rriK 
1 



■x 



(73) 



/2I3 ^ 2 (t2 " tl} + 24^ ~ 3^5 ^ 



d/c4 



-Wfc- 



m. 



2 ^k 



Z^ L-2 ' 



k±0 k2 L 3m| 3 6w k 



E 



J_ [ 2U; k ^k ^7T 

k 2 \-3m n 3 2 



m: 



12oa 



+ 



m: 



12^ 



3 -2(wk-m 7r )(i 2 -ti) 



k^O" " " "* k" 

We distinguished mj{ from f% in the time exponentials even though we consider only the 
case mx = "%. Most momentum sums are divergent and have to be regularized. This will 
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be discussed in more detail when we get to the final expressions including also contributions 
from diagrams (a), (c) and (d). Here we only note that the term proportional to f -^ 
cancels against similar terms from diagrams (c) and (d). 

Of particular interest is the contribution from values of t s between t\ and £2 an d with 
k = 0, which leads to the term proportional to £2 ~t\- Physically, this corresponds to the 
subprocess in which the two pions created by O^ti) are at rest and then scatter into the 
final pions, both also at rest, at any time between t\ and £2- This term can be resummed into 
the time exponentials associated with the out-going pions of the tree-level contribution to 
C(£2, £1) in Eq. (62). (Since there are no such pion-rescattering sub-diagrams in diagrams (c) 
and (d) no similar resummation occurs in those cases.) This corresponds to the finite- volume 
correction to the energy of two pions at rest with isospin 1 = 2: 

E 2n = 2m n + ^2 • (74) 

This volume dependence agrees with the general analysis by Liischer [38] of the finite- volume 
energy shifts of two-particle states. Here we only encounter the tree-level energy shift because 
the pion-pion scattering sub-diagram is a tree-level diagram. 

The terms suppressed by a relative factor exp [— 2(co? — m2)(£ 2 — £^)] correspond to con- 
tributions from excited states. In fact, if m^ > 2m- K (a mass-nondegenerate situation), 
Eq. (50) indicates that the "physical" matrix element that conserves energy-momentum is 
among the nonleading terms. However, in a typical lattice computation, mj{ < 2171^, and the 
finite volume is exploited to obtain the leading term in the large-time limits (£2 ^> £1 ^> 0) in 
which excited states are suppressed. Then, only the unphysical matrix element is obtained. 
We see that one-loop ChPT can be used in order to estimate the contamination from excited 
states, as well as the difference between the physical and unphysical matrix elements (see 
next section). See Appendix B for a discussion of the extraction of the 1 = 2 pion scattering 
length from nonleading time exponentials for diagram (b) in infinite volume. 

25 



The calculation of diagrams (c) and (d) is very similar, and we will not discuss them 
here. In these cases, there are no pion-rescattering sub-diagrams. 

C. Calculation in QChPT 

We now turn to the quenched calculation. As already mentioned in the second section, 
QChPT has an unphysical particle content. Therefore, there is no reason to believe that a 
well-defined Hamiltonian exists, and none of the manipulations of Eq. (50) are justified when 
applied to QChPT. The large-time behavior of Eq. (51) may not be realized. However, a 
Euclidean path-integral type calculation of C(t2, t\) to one loop in QChPT similar to the one 
in unquenched ChPT can be carried out. In particular, the tree- level contribution and the 
term proportional to £2 ~~ ^1 that occurs in diagram (b) are identical to those in full ChPT. 
Resummation of the latter again gives the time exponentials for external legs as shown in 
Eq. (51) with E^ as in Eq. (74). We will extract the matrix element (n^n 10^(0) \K + ) 
from the QChPT result for C{t2,t\) in the same way as in the unquenched case. Note that 
this prescription coincides with the method used in numerical work [8,9,1]. 

One of the differences between ChPT and QChPT is the explicit inclusion of the 77 
field in QChPT. The 77 shows up in QChPT in two ways [11,25]: through arbitrary scalar 
functions of the 77 field multiplying all terms in the Lagrangian and the weak operators, and 
as the ninth component of the nonet field, cf. Eq. (14). It is easy to check that in the case 
of interest, the scalar functions do not lead to new contributions. Furthermore, since the 
operator that represents the chiral currents is 

(L%) it = (SV>„4) . . = ((e»'VAf, S ) 9fl (e-VMf, E t))^ 

*r E9 * Et ),' (75) 

77 couplings would in general arise from diagonal current elements. However, since 

O4 = 2 a 27 [ L 'fj,23 L 'U + ^11-^23 ~ L /j,23 L 22 ~ L ^i22 L 23 + L /j,13 L 21 + L fj,21 L U J ' ( 76 ) 
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the 7/ decouples as a result of the cancellations that happen pairwise between the first and 
the third and the second and the fourth terms (for m u = m^), and the particle content on 
the loops is the same as in ChPT. Since in the mass-degenerate case the double pole only 
occurs in 77 two-point function, no contributions from the double pole show up here. As a 
result, the correlation function C(t2, t\) does not have S and a dependence at one loop. Bad 
chiral behavior caused by double poles is absent. 



D. Collected results from ChPT and QChPT 

The one-loop corrected C(£2,£i) consists of the tree-level term and the one-loop con- 
tributions of diagrams (a) to (d) (for wavefunction renormalizations, see below). Dropping 
contributions from excited states, we obtain for the full and quenched theories, in that order: 



cf(t 2 ,h) 



24ia 27 mlL 3 e - ^* 2- * 1 ) - ™^* 1 



v^/ a 



8ml 



■x 



n m K 



1 + 



m% 



V 



^ + C(^ i ) +B ( m ,L)] + ^ 5 P)(77) 



and 



C 9 (£ 2 ,£i) 



24ialmlL 3 e - ^* 2- * 1 )-"**'* 1 



l 2T lb TT 

V2fi 



Irn^mx 



■x 



mz 



3/2 L8(m^L)3 



C(m n L) + D(m n L) 



in 



(4ttV 2 



hoPW 



where 



C(m n L) 
D{m 1T L) 



{m^L) 



k^O 



7^Er2& + m - 



u k 



2u k ) 
3 



+ 



lW^k_-13^ + 3_^ 
_M3 2^l™_ ,,„ + A, .,3 j ' 



(79) 



[m^L) 6 ^\v(It, 



^k -' ^ k 

and u k is defined in Eq. (72). The contact-term coefficients P and P come from 0(p ) 
weak operators for the full and quenched theories, respectively. Since this result does not 

27 



correspond to the physical weak matrix element, there is no simple relation between P and 
G\ and G<i in Eq. (43). Again, we distinguish between m^ and m^ in the prefactor for 
illustrative purpose only. The term proportional to ti — t\ in Eq. (73) has been resummed 
using Eq. (74). 

In order to extract the unphysical matrix element from C(t2,ti), we need the four- 
pion correlation function Eq. (52) to lowest nontrivial order in ChPT (with / — > f q in the 
quenched case): 

(0|7r + (* 2 )7r°(t 2 )7r-(*i)7r°(ti)|0> = - *** X ~ L% ( l + ^L * 3 ) ■ ( 80 ) 

Ami \ 12r(m„LyJ 

The second term in parentheses is the lowest nontrivial order finite- volume correction for the 

two-pion wavefunction-renormalization. (Note that such terms also occur in Eq. (77) and 

Eq. (78)). We also need the free kaon propagator (Eq. (53)): 

(0|X+(*i)K-(0)|0> = i^^ mKtl ■ (81) 

A cutoff is assumed to regulate the divergent momentum sums in the functions C, D 
and D defined in Eq. (79). The quartic divergences present in the individual terms cancel 
in the combinations C + D and C + D in Eq. (77) and Eq. (78), respectively. The quadratic 
divergence is independent of m n and can be absorbed into a 2 7 or c%. Note that the quadratic 
divergences for the full and quenched theories are not the same. This is another indication 
that the coefficients a 2 7 and a% of the ChPT and QChPT weak operators are different. The 
remaining logarithmic cutoff dependence can be absorbed into the coefficients P and P. 

For large L, we can estimate the size of m^ dependence of the functions D and D by 
replacing the sum j^ X^k °y an integral J ry-ru, with corrections vanishing faster than 
any power of L . The summand in the function C is singular for k = 0, and powers 
of L appear when we convert the sum to an integral. See Appendix C for details. In 
a finite volume, similar momentum sums also appear in the one-loop results for the weak 
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decay constants and wavefunction renormalizations quoted in section 2. These sums are not 
singular near k = 0, and we can replace the sums with integrals as for the functions D 
and D above, with corrections vanishing faster than any power of L . This leads to the 
expressions given in section 2. 

When the appropriate division is carried out in order to extract the unphysical weak 
matrix element, the finite- volume two-pion wavefunction-renormalization corrections (terms 
oc 1/L ) in Eq. (77) or Eq. (78) and Eq. (80) cancel. The matrix element thus obtained is 



renormalized by multiplying it by the wavefunction-renormalization factor Z^yJ r Zx (Eq. (8) 
and Eq. (10) for the full theory and Eq. (19) and Eq. (21) for the quenched theory). We 
also make use of state normalizations Eq. (54) and Eq. (55). The mass-degenerate (m u = 
m^ = m s ) one-loop corrected renormalized unphysical weak matrix elements for the full and 
quenched theories in a finite volume L then are, respectively: 



(<K + TT U \O 4 {0)\K + y = 

24ia 27 m^L 3 / m 2 _ , m 2 



v/2/ 3 V (47T/)* 



61og^f + F( m7T L) + P-K 1 -K 2 



(82) 



and 



[TT + TT \O A (0)\K + ) q 



2 iia^m^L? ( m 2 



v/2/l V (4vr/^ 2 



2 



-31og-f + F( m7r L) + P - K 



up to corrections vanishing faster than any power of L , where the function 



(83) 



17.827 12tt 2 



F{m^L) = —^ + -3 (84) 



is the finite-volume correction calculated in Appendix C. Note that the same finite-volume 
correction occurs in both full and quenched theories. (This is because diagram (b) is the same 
in the full and quenched theories, and because only diagram (b) gives rise to finite-volume 
corrections.) m^ is the measured meson mass. Contact terms with coefficients K\ 2 and K 
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come from wavefunction renormalizations. By comparing Eq. (82) and Eq. (83), aside from a 
change from a 2 -j to a% and / to f q , we see a substantial change to the nonanalytic one-loop 
corrections due to quenching. If we would trade the parameter / for f n in the tree-level 
term of the full theory result, the numerical factor multiplying the chiral logarithm would 
be —15 instead of —6 in Eq. (82), while there would be no change for the quenched result 
since f q = f n up to 0(p ) contact terms in that case. Because of translation invariance, the 
matrix elements of O^fei, 0) are trivially obtained by omitting the factor L from Eq. (82) 
and Eq. (83). 

We end this section with a few remarks about one-loop contributions in the mass non- 
degenerate case (m u = m^ ^ m s ). Due to the mixing of r\ and i], the rj two-point function 
acquires a double pole and hence a and 5 dependence (m u = m^ prevents mixing of it with 
r) and 77). The 77 that appears on the loop of tadpole diagram (a) contributes in a way anal- 
ogous to the way r/-tadpoles contribute to quenched fj{. Hence, we expect the (^-dependent 
term to be singular in the chiral limit, like in Eq. (20). The contribution from diagram (b) 
remains the same as in the mass-degenerate case, because only pions appear on the loop. 
We therefore expect that no "enhanced finite- volume corrections" such as in the case of pion 
scattering [30] at one loop. As for diagrams (c) and (d), virtual mesons of different masses 
appear on the loop. In particular, 77 appears and gives rise to double poles, and hence a and 
S dependence, but no enhanced finite-volume corrections. 

7. Numerical examples and discussion 

Little or nothing is known about the values of the low-energy constants i^, Fj, Gj, P 
and P in Eqs. (30,32,43,82,83). (For some phenomenological information on the F's and G"s, 
see refs. [18,19]; for the others no information is available at all.) Because of this situation, 
we will choose reasonable values of the cutoff A and ignore the contact terms in an attempt 
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to estimate the size of one-loop corrections. Values of the cutoff we will use in the following 
are 770 MeV (the mass of the p) and 1 GeV. We note that in refs. [25,18,19] the value of 
the r] mass was used as cutoff in a phenomenological analysis of 0(p ) contact terms. This 
makes the loop diagrams artificially small, which therefore would most likely underestimate 
the size of the 0(p ) corrections when we set the contact terms equal to zero. 

The determination of the relative size of tree-level and one-loop contributions to the kaon 
matrix elements is particularly sensitive to the definition of the low-energy parameter a 2 -r, 
especially in any scheme that has nonvanishing quadratic divergences. Such divergences are 
subtracted by renormalizing a 2 7, and this introduces an arbitrariness into the determination 
of the relative size of tree-level and one-loop contributions. We will ignore this ambiguity 
in some of what follows, and we will consider a quantity which is independent of a 2 7 at the 
end of this section. 

Apart from the fact that the rj is light in QChPT, the same considerations concerning the 
choice of the cutoff values that apply in ChPT are also applicable in QChPT. Note that the 
cutoff in the quenched theory does not have to be same as that in the full theory. However, 
since non-Goldstone hadron masses in quenched QCD are roughly the same as in full QCD, 
we will again take the cutoff A q for QChPT to be 770 MeV or 1 GeV. In the following 
numerical examples, we have used f^ = 132 MeV, m^ = 136 MeV, mx = 496 MeV, unless 
otherwise stated, m^ has been determined from Eq. (6). 

A. K + — »■ 7r + 7r decay in the real world 

For m n = 0, from Eq. (48), ignoring 0(p ) contact terms, at A = 1 GeV and A = 
770 MeV we find 0.33 + 0.28i and 0.12 + 0.28z respectively for the relative size of the one- 
loop correction to the tree- level term of the K + — ► 7r + 7r matrix element at m^ = 0. This 
is smaller than what one would obtain from the results of ref. [16], which (for the real 
parts) would give 0.57 and 0.35, respectively. This is due to the fact that in Eq. (48) we 
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kept contributions from the nonanalytic terms that become numerical constants after setting 
run = 0. This example demonstrates that estimates of the relative size of one-loop corrections 
are quite sensitive to the choice of the values of 0(p ) low-energy constants, or, equivalently, 
the cutoff, for quantities involving kaons. In addition, this matrix element is quite sensitive 
to the value of the pion mass, as one finds by substituting m^ = 136 MeV in Eq. (43) with 
0(p ) contact terms ignored, which gives considerably larger relative one-loop contributions 
for the real parts: 0.63 + 0.20i and 0.36 + 0.20i respectively. (Br , in Eq. (36), is much less 
sensitive to the value of m^.) 

It is instructive to estimate the change in magnitude of 0(p ) contact terms due to 
a change in cutoff. This change can be calculated from Eq. (43), given two values of the 
cutoff. For the m^ = case, a change of A from 1 GeV to 770 MeV induces a change 
AGi = 2.35. For the m^ = 136 MeV case, one finds A(G\ + 2/7^2) = 2.96. Here, we defined 
G\ — G\ + -§-K\ — j*K2 and G2 = G2 — ^K\ + ^K 2 - Since the changes in G\^ 2 are of order 
one, we believe that using values for the cutoff of 770 MeV and 1 GeV gives a reasonable 
estimate of the uncertainty in the size of one-loop corrections coming from our ignorance of 
the values of the 0(p ) low-energy constants [39]. 

B. K + — ► 7r + 7r decay on the lattice 

The physical and unphysical K + — »■ tt + tt matrix elements in the full theory can be 
related up to one loop using Eq. (43) and Eq. (82): 

(, + Aom\K + ) phys = x ^|-^ (- + Aom\K + ){ nphys , &) 

IT 

where 
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x = Stjt^t m : t^ (86) 



where the subscript "(un)phys" refers to the (un)physical matrix element (we imagine the 
unphysical matrix element to be obtained from a lattice computation), and the superscript 
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"/" refers to the full theory. Eq. (85) relates the real-world matrix element to the mass- 
degenerate unphysical matrix element as it would be computed on the lattice with degenerate 
quark masses. The relation is obtained by reconverting the factor l//£ back to 1// in 
Eq. (43). The one- loop corrections in the resulting expression after this change are those 
given in Eq. (43), but without If, and now the 0(p ) contact terms are as in the numerator 
of Eq. (86). The parameter U, then, is the ratio of the nonanalytic one-loop contribution 
and the tree-level contribution in Eq. (43), leaving out If. 



m K T m K 



1 ' = 7TTV \ h + 2 %K 2 ( ] " + l >> + [ <-+<i ] \ ~ ( s7) 

\4*U) \ m z K -m% 

and equals 0.089 or —0.015 for A = 1 GeV or A = 770 MeV respectively (we will ignore the 
imaginary part of U, since it does not contribute to the magnitude of the amplitude to order 
p ). We use upper-case symbols F n and M n for the decay constant and meson mass as they 
would be computed on the lattice, and lower-case symbols for the real-world values of these 
quantities. At tree level, X — 1. Note that the ratio X is independent of a 2 7, so that it 
is not sensitive to ambiguities in the renormalization of this parameter. Note also that the 
0(p ) contact terms do not cancel between the numerator and denominator of Eq. (86). 

If we take values for the pion mass and the volume typical for lattice computations, 
for example M w = 0.4 GeV, M n L = 8, and F w = f n (and again set the 0(p ) low-energy 
constants equal to zero), we find the denominator of X to be about 1.8 (leading toX« 0.6) 
for A = 1 GeV, and about 1.6 (again leading toXsi 0.6) for A = 770 MeV. 

Actual lattice computations have all been done in the quenched approximation [8,9,1], 
and it is therefore interesting to compare the full physical and quenched unphysical (mass- 
degenerate) matrix elements directly. In this case, we find, using Eq. (43) and Eq. (83) 



, |O4(0)|A' + )„, 1 „ s = Y ^ (41 ! %7^ (* + AO 4 (0)\K + )'! 
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(89) 

9 1^^ 1vj -tt 1 ZTV/l/T T\ 1 D ^ 

where A q is the cutoff in QChPT. Again, at tree level Y — 1. In this case there is no reason 
that a 2 7 = «27 ( or / — fq)i since we are not comparing matrix elements in the same theory. 
In particular, Eq. (88) is sensitive to renormalization ambiguities in these parameters. 



Lattice parameters 


Ml 


M n L 


^(1) 


v (0.77) 

1 (0.77) 


1 (0.77) 


v (0.77) 

^(1) 


16 3 x 25(or x 33) 

P = 5.7 or 

1/a = 1.0 GeV 


0.168 
0.263 

0.474 


6.55 
8.20 
11.0 


0.72 
0.68 
0.65 


0.69 
0.67 
0.70 


0.77 
0.75 
0.77 


0.65 
0.61 
0.59 


24 3 x 40 

f3 = 6 or 

1/a = 1.7 GeV 


0.303 

0.428 


7.77 
9.24 


0.65 
0.63 


0.66 

0.67 


0.72 
0.74 


0.59 

0.57 



Tafr/e 1. T/ie factor Y, Eq. (89), with contact terms ignored, at different values of M^ for 
different combinations of values of the cutoff A and A q for the full and quenched theories 
respectively. Lattice parameters including the four-volume L xT are shown. The superscript 
on Y denotes A in GeV; the subscript on Y denotes A q in GeV. M* is in GeV . 

The tree-level relation (Y = 1) was used in refs. [9,1,4] in order to extract the physical 
weak matrix element from the unphysical one computed on the lattice, assuming a 2 7 = ^27 
and f = fq. It was found that the lattice values, corrected in this way, were at least 50% 
larger than the experimental value. It is therefore interesting to estimate the values of the 
one-loop correction factor Y. We will again set the 0(p ) low-energy constants equal to zero, 
and consider two values of the cutoff, 1 GeV and 770 MeV, for both the full and quenched 
cases. Since these are different theories, there is no reason that the cutoff should be chosen 
the same in the numerator and denominator of Eq. (89), and we will compute four different 
values of Y for each M n corresponding to the four possible different choices of the cutoffs in 
both theories. We will take the spread in the values of Y obtained in this way as an estimate 
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of the uncertainty in Y . Using F n = f n , and values for M^ and L from refs. [4,1], we report 
the values of Y summarized in table 1. 

In order to get an idea about the size of finite- volume effects on Y, we list in table 2 the 
values of Y without the finite- volume correction F^M^L) in Eq. (89). 



Lattice parameters 


Ml 


^(1) 


(0.77) 

(0.77) 


y(l) 
(0.77) 


(0.77) 

r (l) 


16 3 x 25(or x 33) 

(3 = 5.7 or 

1/a = 1.0 GeV 


0.168 
0.263 

0.474 


0.82 
0.79 
0.79 


0.80 
0.80 
0.88 


0.88 
0.88 
0.98 


0.74 
0.71 
0.71 


24 3 x 40 

P = 6 or 

1/a = 1.7 GeV 


0.303 
0.428 


0.78 
0.78 


0.81 
0.86 


0.89 
0.95 


0.71 
0.71 



Table 2. The factor Y , Eq. (89), with both contact terms and finite-volume correction 
F(MkL) ignored. For further explanation, see the caption of table 1. 

It is evident from a comparison of tables 1 and 2 that the finite-volume correction can 
be substantial, of order 10-20 %. Note that the finite-volume corrections do not take into 
account any corrections vanishing faster than any power of L . While these are very small 
at large values of M^L, they may be substantial for the values considered here [30]. We 
do not address this issue further in this paper. Of course, it should be stressed that these 
estimates of the correction factor Y do not take into account the factors a^/oli-j and (fq/f) 
in Eq. (88). The fact that quadratic divergences differ between the full and quenched theories 
(cf. section 6.D) may be a hint that the ratio c^/a 27 is n °t very close to one. 

The values of Y in tables 1 and 2 were evaluated by first calculating the numerator and 
denominator of Eq. (89), and then dividing the results. We could also have first calculated 
the 0(p ) terms in the numerator and denominator, and then evaluated Y as one plus the 
difference of these 0(p ) terms, since, in order to distinguish between these two methods, 
a two-loop calculation would be necessary. Since the values for Y in tables 1 and 2 differ 
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from one by rather substantial amounts, we would obtain significantly different results for 
Y . In fact, the average values of Y would have been lowered by about 10-20 % (relative to 
one). This difference reflects the uncertainty introduced by cutting off the chiral expansion 
at one-loop order. 




Figure 2. The physical K + — ► 7r + 7r° amplitude vs. M% from tree-level and one-loop 
(Q)ChPT. The one-loop corrected points lie below the tree-level corrected points. For de- 
tailed explanation, see text. Squares and octagons refer to data with four-volume 16 x 
25 (or x 33) and (5 — 5.7 and four-volume 24 x 40 and (3 = 6 in table 1, respectively, 
c = 2\pl /(Gpsin6 c cos6 c ). 



In fig. 2, we replot some of the lattice data in fig. 16 of ref. [4] (or fig. 3 of ref. [1]) 
for the K + — ► n + n matrix element at different values of M% (errors are statistical). We 
chose data with M n < 770 MeV and with the larger spatial volume for f3 — 6. The M^ 
and L values are those used in tables 1 and 2. These points have already been corrected 
for the difference between physical and unphysical values using tree-level ChPT, assuming 
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«27 — a 27 an d f — fq- Again assuming these equalities, we also plot the one-loop corrected 
physical matrix element using Eq. (88) with values of Y as shown in table 1. In all cases, 
the one-loop corrected points lie below the tree-level corrected points. At each point, the 
central value is obtained using the average of the four different values of Y, while the error 
bar is obtained from the largest and smallest values of Y. These error bars do not include 
the statistical errors of the lattice computation. We see that the one-loop corrected values 
are closer to the experimental value. 

We end this sub-section with some discussion of the 0(p ) low-energy constants and the 
contributions from excited states to the correlation function Eq. (49). 

First, we note that, for the quenched unphysical weak matrix element Eq. (83), a change 
of cutoff A q from 1 GeV to 770 MeV results in a shift A(P — K) = 1.57, again of order one. 

Due to the fact that quenched and full QCD are different theories, there is no relation 
between the low-energy constants P, K and G\ 2, K\ 2, and therefore, in principle, no 
cancellation occurs in Eq. (89). Alternatively, at small enough M n , one may envision relating 
physical and unphysical weak matrix elements, both in the full theory, as in Eq. (85) when 
realistic lattice computations for full QCD become possible. However, since the unphysical 
matrix element is calculated at different external four-momenta, the contact terms P for the 
unphysical full matrix element, Eq. (82), and G\ 2 again have no simple relation, and again 
no cancellation takes place in Eq. (86). While in principle the constants P, K, P, K\ and 
Ki can be obtained from the lattice, this is not the case for G\ and G^ as long as only 
matrix elements with all mesons at rest are computed numerically. The values of coefficients 
K\ and K^ associated with weak decay constants have been constrained in ref. [25] (where 
a different notation is used) with the help of large- N c arguments. A larger number of 0(p ) 
operators contribute to K + decay, and in Eq. (43), G\ and Gi are actually particular linear 
combinations of the associated coefficients. A determination of some of the coefficients in 
the m^ = approximation has been attempted in refs. [18,19]. 
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A determination of the prefactor (fq/f) requires a more precise knowledge of weak decay 
constants both in full and quenched theories in the chiral limit than appears to be available 
at the moment. A fit of quenched lattice data for F n in the mass-degenerate case to Eq. (18) 
would give f q and K [40], but /, K\ and K<i are only poorly known from phenomenology 
[25]. In addition, nothing is known about the ratio 0:27/0:27 • As long as no more detailed 
knowledge of all low-energy constants involved is available, we are limited to the estimates 
of ChPT one-loop effects as presented in this paper. 

In order to separate the contribution of the lowest two-pion state (for which each pion 
has spatial momentum k = 0) from that of the closest excited states (for which each pion 
has spatial momentum |k| = 2ir/L (cf. Eq. (73))) to Eq. (49), we demand that 



a = 2M^l + (^) 2 -l)»l. (90) 

where t = t% — t\ corresponds the lower bound of the fitting range used in a numerical 
computation. For lattice data listed in table 1, fitting ranges in the four- volumes 16 x 
25(or x 33) and 24 x 40 are about 5 ~ 9 and 9 ~ 13 in lattice units, respectively [40]. As 
an example, for the smallest mass M£ = 0.168 in table 1 we find f2 = 1.58 for t = 5, and 
condition Eq. (90) is not really satisfied. For all mass values in table 1, Q ranges from 1.04 
to 1.67. Note that in general, at sufficiently large M n L and fixed L and t, Q decreases as 
M n increases. Since there are six excited states with |k| = 2ii/L, larger values of Q would 
be needed to be confident that the excited states are not contaminating the results. 

C. Relating (7r+7r |O 4 |K + ) and (K°\O'\K ) 

An interesting ratio is that of (tv^tt \0/±\K + ) and (K \0'\K ), because it is independent 
of a 2 7 (or a 27 in the quenched case). We will consider the dimensionless quantity 

<^W> (91) 

(k u \o>\k0) 
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both for the physical and the unphysical case, where in the latter case the unphysical weak 
matrix element (ir + % \0^\K + ) is used. This ratio was examined in ref. [21] at tree-level for 
the physical case. We can parametrize 1Z as 

f Si 171%- — ml , N . , 



for the physical case, and 



unphys /9 \ unphysj \ ' 



for the unphysical case (with degenerate quark masses). In these equations, R stands for 
one-loop contributions, i.e. R = at tree level. (Note that, at tree level there is no difference 
between the full and quenched theories in the unphysical case.) 

At one loop, we evaluate R p hys again at two choices of A: 1 GeV and 770 MeV, setting 
0(p ) contact terms equal to zero. Using Eqs. (7,9,30,43) we find that R p h ys equals —0.94 + 
0.20i and — 0.62 + 0. 20i at A = 1 GeV and 770 MeV respectively. These are large corrections, 
and call into question the reliability of ChPT for this case. 

In the unphysical (mass-degenerate) case, we find, from Eqs. (7,9,30,82), 



1 / M l 2 \ 

^2 (31og-^ + F(M W L) + P - F l + -(K 1 + K 2 )j 



unphys _ ^ ^f __ ^ M n („, M^ 

3i/y/2 unphys {^F^ 

(94) 

for the full theory, and, from Eqs. (18,20,32,83), 

5^ = i + jR 9 =1 + ^rf( 3log Ml + F (M 7r L)+P-F 1 + -K) (95) 

3i/y/2 unphys (4nF n ) 2 I & A 2 q K J i 3 y / v; 

for the quenched theory. Note that the nonanalytic terms are the same in Eq. (94) and 
Eq. (95). 

For numerical estimates of the one-loop corrections, we ignore the contact terms and set 
F-k = f-Ki as before. Table 3 lists the resulting one- loop corrections for the same lattice data as 
in table 1. Both one- loop corrections in the infinite- volume limit (the chiral logarithm only), 
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labeled as R°°, and with finite-volume corrections (F(M 7T L)), labeled as R , at different 
cutoffs are shown. 



Lattice parameters 


Mi 


poo 

R (l) 


R (l) 


poo 

-""(0.77) 


T3L 

-""(0.77) 


16 3 x 25(or x 33) 

(3 = 5.7 or 

1/a = 1.0 GeV 


0.168 
0.263 

0.474 


-0.33 
-0.38 
-0.39 


-0.14 
-0.15 
-0.092 


-0.23 
-0.23 
-0.12 


-0.039 

-0.0050 

0.18 


24 3 x 40 

(3 = 6 or 

1/a = 1.7 GeV 


0.303 
0.428 


-0.39 
-0.40 


-0.11 
-0.072 


-0.22 
-0.15 


0.058 
0.17 



Table 3. Nonanalytic one-loop corrections in Eq. (95) in infinite volume (R°°) and in finite 
volume (R ) are shown. The subscript on R denotes A q in GeV. M^ is in GeV . 

Large positive finite-volume corrections partially offset the negative corrections from 
chiral logarithms. The uncertainties for the ratio Tt , associated with the range of 

° unpnys ° 

cutoff explored are roughly 10% to 15% for the lower two values of M n , and increase with 

8. Conclusion 



In this paper we have presented an analysis of the nonleptonic kaon decay mode K + — ► 
7r + 7r to one loop in chiral perturbation theory, with emphasis on the comparison of results 
obtained from lattice QCD computations with experiment. One-loop ChPT makes it possible 
to gain analytic insight into three different systematic errors: quenching effects, finite-volume 
effects, and effects from the fact that lattice computations (so far) are at degenerate quark 
masses and have all external momenta vanishing. The last effect implies that on the lattice 
the relevant weak matrix element has been computed at an "unphysical" point. We have 
found that, for the most recent lattice results [1], all three effects are important. 
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Refs. [1,4] already employed tree-level ChPT in order to compare lattice results with 
experiment. At tree-level no quenching or finite-volume effects show up in ChPT, while 
the conversion factor between the unphysical lattice matrix element and the physical one 
amounts to a simple mass ratio (Eq. (85) with X = 1, or Eq. (88) with Y = 1, / = f q and 

a 27 = aj 7 ). 

This tree-level conversion factor between unphysical and physical matrix elements gets 
corrections at one loop, which are embodied in the one- loop correction factor X in Eq. (85). 
X is given by Eq. (86), which also includes finite- volume corrections. We found that, for 
a meson mass M^ = 0.4 GeV and a volume M^L = 8, X is approximately 0.6, i.e. a 
substantial correction. (At M n L = oo, one finds X rs 0.7.) 

Eq. (88) gives the conversion factor between the unphysical, quenched matrix element in 
finite volume and the real- world, physical matrix element. The one-loop effects are contained 
in the factor Y, given in Eq. (89). Table 1 gives estimates of the factor Y for meson masses 
and spatial volumes used in ref. [1], while table 2 gives Y for the same masses, but in infinite 
volume. From a comparison of these two tables, we see that finite-volume effects cannot be 
ignored for these masses and volumes. In order to isolate the effect from quenching on Y 
we can take the ratio of Y and X. Considering only the chiral logarithms in this ratio, one 
finds that Y/X « 1.2 for M n = 0.4 GeV, F n = f n and A = A q = 1 GeV (finite-volume 
corrections cancel between Y and X). 

In refs. [1,4], lattice results for the K + — > 7t + n matrix element were corrected using 
tree-level ChPT, and compared with experiment. Some of the data points of ref. [1,4] are 
reproduced in fig. 2. In fig. 2 we also showed how the lattice points would shift, if our best 
estimates for the one-loop correction factor Y are included (for the error bars, see section 
7). From fig. 2, we conclude that one- loop effects in ChPT narrow the gap between lattice 
results and the experimental value. 

We should stress that there are other possible sources of discrepancy between lattice 
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results and the experimental value. First, in obtaining the one-loop corrected points in fig. 
2, the factors a^ifdi-j and (fq/f) in Eq. (88) have been ignored. Since the quenched theory 
is different from the full theory, there is no a priori reason that these factors should be equal 
to one. Furthermore, our estimates of the factor Y have been obtained by choosing some 
values of the cutoff in (Q)ChPT, and ignoring all 0(p ) low-energy constants. While the 
error bars on our results in fig. 2 give some indication of the uncertainty introduced by 
this procedure, more information on these low-energy constants would be needed in order 
to be more confident about our estimates of Y . Third, since one-loop effects turn out to be 
rather substantial, one may expect that two-loop corrections cannot be completely neglected. 
Also, our estimate of finite-volume effects included only effects that go like inverse powers 
of L, and do not take into account exponential corrections. (For more detailed discussion of 
some of these points, see section 7.) Finally, there are possible systematic errors associated 
with approximations that we did not discuss at all, such as contamination of lattice results 
by higher excited states in Eq. (49), breaking of isospin symmetry and scaling violations. 
However, we believe that we have demonstrated that one-loop corrections in ChPT do move 
the lattice results closer to the experimental value. 

If quenched lattice computations were performed with nondegenerate quark masses, this 
would introduce new systematic effects into the difference between lattice results and the 
experimental value of the K + — *■ 7r + 7r matrix element, which would show up at one loop 
in QChPT. With nondegenerate masses, the quenched result would be dependent on the 
"■//-parameters" S and a, and presumably it would diverge in the chiral limit (cf. section 
6). Experience with other quantities in QChPT suggests that these new effects could be 
rather substantial, which would imply that nondegenerate quenched computations may be 
of limited value. 

In section 4, we presented our one-loop calculation of the unquenched, physical K + — ► 
7r + 7r matrix element with nondegenerate masses. For m^ = 0, our result agrees with that 
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of ref. [16]. However, we note that this quantity is particularly sensitive to the pion mass, 
so that its numerical value at the physical pion mass can deviate significantly from that at 
ititt = 0. 

Finally, we reviewed and extended previous existing one-loop calculations in ChPT of 
the kaon B-parameter, B^ [16,36,19,29,20,5] (cf. section 3). We extended the quenched 
calculation of refs. [29,5] to include the dependence on the 77 -parameter a. Furthermore, we 
discussed a modified version that could be of use for staggered fermion QCD (cf. Appendix 
A). In section 7, we used our results in order to calculate a quantity which is basically the 
ratio of the K — K and K + — > 7r + 7r matrix elements. This ratio is independent of the 
0(p ) low-energy constant a 2 7, and therefore is of some interest (it is not independent of 
0(p ) low-energy constants). For the real world, we found that this ratio has large one-loop 
corrections (mostly due to large chiral logarithms in the K — K matrix element). For the 
lattice data of ref. [1] these one-loop corrections change the tree-level value by some +20 to 
-40 % (cf. table 3). 
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Appendix A: B'L with staggered fermions 

In this Appendix, we present a somewhat modified expression for B'L which may be 
relevant for quenched staggered-fermion lattice computations of this quantity at values of 
the lattice spacing where the continuum flavor symmetry is not yet completely restored. 
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It has been observed [41,42] that the double-pole contribution from the Goldstone mesons 
appearing on the loops in QChPT have to be (staggered-) flavor singlets, while the external 
mesons are usually taken to be the one exact Goldstone boson (per staggered fermion field) 
of staggered-fermion QCD, which is not a flavor singlet. At typical values of the lattice 
spacing, one observes that the nonexact Goldstone bosons are more or less degenerate [43], 
while their masses do not go to zero for zero quark mass, due to the staggered breaking 
of chiral symmetry. Of course, (Q)ChPT, being a continuum effective theory, is inherently 
inadequate to accommodate scaling- violation effects. However, heuristically, if we denote the 
masses of the nonexact Goldstone bosons as fnx an d tu-k-, B'L will depend on these nonexact 
Goldstone masses through the double-pole contribution to the one-loop corrections. Inserting 
mj{ and m^ into the double-pole terms, and ignoring 0(p ) contact terms, we now get 
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(96) 



2e ° 1 - e 

with e = 1 — m^/mj^. In this expression, mj{ and e refer to the exact Goldstone boson 
masses. Note the appearance of terms that go like l/m K . We should stress that this 
expression is heuristic, since it assumes that the only effect of nonzero lattice spacing is the 
appearance of nonexact Goldstone masses in the double-pole terms. 



Appendix B: Discussion of nonleading 
time exponentials from diagram (b) 



In ref. [22] Maiani and Testa showed that, for any decay process with two identical final 
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pions as outgoing particles, the leading exponential time dependence for the final pions in 
the corresponding three-point Euclidean correlation function acquires power corrections in 
the special case that both pions are at rest, and when first the infinite- volume limit is taken. 
The leading correction comes from the sum over all excited states in Eq. (73) (cf. last line 
in that equation), which for infinite spatial volume becomes an integral. It has the form (cf. 
Eq.(lO) ofref. [22]) 



e -W»-'i) -> e -W*2-*i) / x _ fl J=2 / ™* J . (97) 

Note that the second term in Eq. (97) differs by a factor 1/2 as compared to Eq. (10) in ref. 
[22]. This takes into account the phase-space difference between distinct and identical pions 
in the final state. The parameter Qq is the 1 = 2 pion-pion ,S-wave scattering length. Note 
also that E^ —* 2m T: in the exponent, since this is an infinite- volume result. Indeed, if we 
replace the sum over k in the last line of Eq. (73) by an integral, we find Eq. (97) with 



1=2 _ 1 m -n 



4 =-—-£< (98) 



which is the tree-level ChPT expression for the 1 = 2 S'-wave scattering length. This infinite- 
volume result constitutes a useful check of our calculation. 

Appendix C: Finite-volume expansion of C(m n L) in Eq. (79) 

In ref. [38] Liischer proved that, for a function / of which all derivatives are square- 
integrable, 

with integer q > and f^'(0) = ( -jH ) f{k )|^.2_q, up to corrections vanishing faster than 
any power of L . A^ is the Laplacian with respect to k and z(q) = Zoo(q,0) is a zeta 
function [38], with 

z (\) = -8.91363292 , z(0) = -1. (100) 
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For a smooth function g(k ) for which the sum in Eq. (99) is ultraviolet divergent, we 
introduce a regulator mass M, and apply Eq. (99) to 

/ M 2 \ n 

f^=^[WTW-) • (101) 

where the positive integer n is chosen sufficiently large to make / square-integrable. For 

M 2 ^oo, 

/(k) — > «?(k) , (A fc )(^ 1 )/(k) — > (A fc )(«" 1 )(/(k) and /%) — > /?)(0) . (102) 

We see that the finite- volume corrections in Eq. (99) are finite in the M — »■ oo limit, 
and could have been calculated directly from g, while the integral "carries" the ultraviolet 
divergence. Applying this to Cim^V) in Eq. (79) and keeping only the universal logarithmic 
cutoff dependence from the integral term, we find 

^, n 3 , m 2 3,2(1) 1 9 1 , ino , 

Cim^L)-*— — -rrlog^ 1 hr r + - : r~ , (103) 

V ; 2(4tt) 2 & A 2 8tt 2 m^L 4 (m^L) 3 ' V ; 

valid up to corrections vanishing faster than any power of L . For a discussion of the 
power-like ultraviolet divergences and contact terms we refer to section 6. 
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